BACKGROUND OF THE INVENTION 

deleterious effect, of that ,^ «d how to combat the 

°SSr ^^ ^^ Um ^ l ° m ada ^ opric 

vw W^2^ fa** each with a kno™ diversity rphasc 

Gowriv* and others to deduct flwh 1 u^m fcind 811(3 the m « hod 
embodim** of this a.thod " ' - Tdesco P s In * e ^ 

Reusing the optical a ^tc ph f e shift which can be induced b v 

v-fach use phase diversity. " ' ^ Ct al " and ?axma, i hold current oaten?* ' 

i-gc scjui^ this invention uses on,y the ^ured 

required. Successive ima^LmthlSiS^ ^ 85 * dcft ^ in«ge is 

and the change, *> the ad^tp^tt^ * 

adapts S T™£^:? C ™* f - .echanfcn*. It applies to any 

tbe diversity in^e* and iffS of T 1 the use of sequential images « 

optic. ' ie 01 a st{ ^ n <^ diversity processor to control the adaptive 

Referenced baCicground * fcreaces « ^ted below 

Plane I^uSs^ " Phase ^ Image and Difixacti on 

J ^"^^ Soc. rfAa.66,961 (1976, 

5. B. Ellerbrock a Mo£ on 1 M ^'* °P tics "' <*• S"" 21. 829 (1982) 

90 (1982J». Mor "^', Linear methods m phase retrieval/ 1 SPIE Proc. 35 1, 

? " AmXt V l6o'aS7) RCtriC " 31 by DifferCntial Measurements", J. Opt. Soc of 

10. IB^eT^v' ASt f r miCaI Sode * ° f dlc Pa * 103, 13 099:). 
c'v^:;.- S^^ f Sff PtiVH ^ m The 





BRIEF DESCRIPTION OF THE DRAWINGS: 

adaDtivfltl 1 '! a . Wmbined .^Pictorial representauon of an imaging system which uses an 
adaptive opto and a sequent^ diversity processor, in accordance with the invention. 

« V g " IS * block of the sequential diversity processor, showing the input 

images from a detector and the output control signals for the adaptive optic. 
. ««■ * m * a con * uter simulation of the invention as it would apply to rhe imaging of 

a point source .like a star, imaged through a turbulent atmosphere. 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

The detector seen in both Fig. 1 and in Fig. 2 is one in which the picture is captured as 
discrete picture elements - essentially, a digital camera. We propose no new detecting 
mechanism fur the camera, digital cameras exist and are improving in sensitivity and 
resolution at a dramatic pace. 

The adapt*** optic (AO) shown in both Fig. 1 and in Fig. 2 is a high-resolution device 
which allow a wide range of correction mechanisms, such as Zernike polynomial firing of a 
complicated wavefront. It is similar to the sophisticated systems currently installed on national 
telescopes and tt has capabilities well beyond the automatic focussing and registration 
mechanisms in consumer camcorders. These sophisticated AO systems improve image quality 
by about a factor often. They will almost certainly appear in consumer products, after the cost 
ot adaptive optics decrease and after simple, real-time control mechanisms are invented. 

^ The AO in these figures is used to cancel the wavefiont distortion, W. Thus, to be 
effective, the optical system of Fig. 1 must sense the wavefiont There are five well-known 
mechanisms for wavefront sensing, which we review here. The five are as follows: 

1 • Dithering. This method continuously changes the adaptive optic and monitors 
the image quality of the observed image. The method was patented by O'Meata, who made 
slowly varying changes to each channel of the adaptive optic, modulated each channel, 
demodulated the image quality in each channel, and controlled the adaptive optic with the set 
of image qualities. A recent improvement, stochastic parallel-gradient-descent, has been 
proposed by Voronstov and Sivokon (Reference 12). They make random changes in all 
channels and, over time, determine which average changes improve the image quality . 

2. Shearing Interferometer. This method uses a reference beam to create an 
interference pattern which describes the unknown wavefront. Hardy was the first to patent a 
wavefront sensor based on this device. It has been used successfully in several astronomical 
installations. Iu, major drawback is the extra measuring equipment, including a laser-based 
interferometer, that has to be installed. 



3. Shack-Hartmaan Sensor. This device uses an array of lenses to image small 
sections of the wavefront onto a detector. Shifts in the small images are caused by local tilts in 
the waveform. From these tilts one reconstructs the wavefront. Feinleib held the first patent. 
A major drawback is the requirement for the lenslet array and sensors. 
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4. Curvature Sensing, This uses two or more images, measured along the path of 
the optical system. The local curvature cf the propagating wave is determined and it is 
propagated, by computer simulation, back to the aperture. The computer-generated phase 
controls an adaptive optic, The method was proposed by Stoddier and related patents were 
awarded to GafTard. et. a! and Rafanelli, et al„ A system of this type is in use on a telescope in 
the state cf Hawaii. References 6, 7, 9 and 10 describe systems of this kind. 

5 . Phase retrieval and diversity. Our invention is similar to these, in that they use 
digital imagery, an adaptive optic, and multiple images. The method was first suggested by 
Gonsalves (Reference 2). References 1-5, 8 and U use phase retrieval and/or diversity. 

Many inventors have tried to use post-processing of a series of images to enhance The 
composite. Patents to Rhode* fall into this category. While such techniques may be useful if 
the data are already captured and no adaptive optic was available during measurement, this 
inventor teaches, rt If ycu want a good picture, take a good picture. 11 This implies that it is far 
better to improve ihc image in its taking than to try to improve the image by computer 
processing, after the fact. 

Referring to Fig. 1 and Fig, 2, we use the following notation: 

Wfk) m Unknown distorting wavefiront at time k, 

T(k) - Phase put on the AO at time k. 

C(k) ~ Compensated phase to be estimated by a diversity algorithm 

- \v(k> +T(k). a: 

I'M) 33 Measured image at time k. 

D(k) = Diversity phase. 

We consider a diversity algorithm where I(k-l) is the first image and I(k) is the 
diversity image. With this convention the diversity phase is the change in the AO phase from 
time k- 1 to time k. Thus, 

D<k) - T(kUT(k-l). (2) 

The phase diversity algorithm is set up to estimate C(k), the compensated phase at time 
k Call the estimate Q(k). Q(k) is, from equation (1), 

Q(k) - Wi(k) + T(k) f (3) 

where Wl(k) is an estimate of W(k) f the unknown pha$e 2t umc k. (We know T(k). so it need 
not be estimated.) At time k ^ 1 we would like to set Iht AO phase to the negative of the 
unknown wavefiom at time k+I . We do not know W at time k+1 but eve do have an estimate 
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we set 

T(k+1) = - Wl(k) , 
which will tend to cancel the wavefront distortion at k+1, 
Solving (3) for Wl(k) and substituting it into (4), we have 

T(k+1) - - Q(k.) +T(kj , 
which implies 

T(k) - -Q(k-1) + T(k-1). 



(5) 



Turning to the diversity phase, we put (5) into (2) to get 

D(k) -( - Q(k-t) + T(k-l)) - T(k-1) 

= -Q(k-0. C 6) 
This is the specification for the diversity phase at time k. 
Finally, we put (6) into (5) to get a new expression for the AO phase: 

T(k) = T(k-1) + D(k) . (7) 

Equations (6) and (7) gi ve us the recipe we sought. Hie AO phase, T(k), is the previous AO 
phase plus the diversity phase, D(k); and the diversity phase is the negative of the diversity 
algorithm's output at time k-1. 

These equations are expressed in block diagram form in Fig. 2. The inputs to the 
diversity algorithm are Images I(k-1) and I(k); the output is Q(k), an estimate of the 
compensated phase, the phase diversity is D(k); and a feedback loop calculates T(k). 

Fig . 3 shows a computer simulation of 24 cycles of the sequential diversity imager. 
The object is a point source, such as a star. The first observed image of the star is in the lower 
left of Fig 3 . The compensated image is directly above. The last pair is in the upper right. 
The lower rows show how the star "twinkles'' over time. Each uncompensated image is called 
a speckle pattern and is characteristic of a star imaged through the earth's atmosphere. 

A measure of the image quality is the Strehl ratio, i. e., the peak value of the aberrated 
star image to the peak value of an unabbeiated star image. Higher Strehl ratio means higher 
image quality. For the 24 images shown, the average Strehl ratio is 0.045 . 

The sequential diversity imaging technique produced the images shown directly above 
each twinkling image. The average Strehl ratio for the compensated images is 0.467. This is 
an improvement of about a factor of 10. 
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